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Abstract 

In a simple holographic model for QCD in which the Chern-Simons term is incorporated to 
take into account the QCD chiral anomaly, we show that baryons arise as stable solitons which 
are the 5D analogs of 4D skyrmions. Contrary to 4D skyrmions and previously considered 
holographic scenarios, these solitons have sizes larger than the inverse cut-off of the model, 
and therefore they are predictable within our effective field theory approach. We perform a 
numerical determination of several static properties of the nucleons and find a satisfactory 
agreement with data. We also calculate the amplitudes of "anomalous" processes induced by 
the Chern-Simons term in the meson sector, such as a; ^ 7r7 and lv — > Stt. A combined fit 
to baryonic and mesonic observables leads to an agreement with experiments within 16%. 



1 Introduction 



In the large- Nc limit, strongly interacting theories such as QCD have a dual description in terms 
of a weakly-interacting theory of mesons [1]. In this dual description, baryons are expected to 
appear as solitons made of mesons fields, usually referred as skyrmions [2,3]. 

Skyrmions have been widely studied in the literature, with some phenomenological successes. 
Nevertheless, since the full theory of QCD mesons is not known, these studies have been carried 
out in truncated low-energy models either incorporating only pions [2,3] or few resonances [4]. 
It is unclear whether these approaches capture the physics needed to fully describe the baryons, 
since the stabilization of the baryon size is very sensitive to resonances around the GeV. In the 
original Skyrme model with only pions, for instance, the inverse skyrmion size equals the chiral 
perturbation theory cut-off A^pt ~ 47rF7r, rendering baryon physics completely incalculable. Other 
examples are models with the p-meson which were shown to have a stable skyrmion solution [5]. 
The inverse size, also in this case, is of order rrip ~ A^p^, which is clearly not far from the mass 
of the next resonances. Including the latter could affect strongly the physics of the skyrmion, or 
even destabilize it. 

In this article we will consider a very simple five-dimensional model for QCD, which has already 
been shown to give a quite accurate description of meson physics [8-10]. This 5D model has a cut- 
off scale A5 which is above the lowest-resonance mass rup. The gap among these two scales, which 
ensures calculability in the meson sector, is related to the number of colors A^^^ of QCD. In the large 
Ac-limit, one has A^/rrip cxd and the 5D model describes a theory of infinite mesonic resonances. 
We will be interested in studying the solitons of this 5D theory that will correspond to the baryons 
of QCD. We will find that these 5D skyrmion-like solitons have an inverse size p^^ ~ nip smaller 
than the cut-off scale A5. Therefore, contrary to the 4D case, they can be consistently studied with 
our 5D effective theory. The expansion parameter which ensures calculability will be provided by 
l/(pA5) < 1. 

Once calculability is established, it is meaningful to compare our predictions with experiments. 
We will study numerically the 5D baryons and calculate several static properties of the nucleons 
such as the axial coupling, magnetic moments and radii. An important ingredient of our model 
will be the Chern-Simons (CS) term that not only will incorporate the QCD anomaly, but will 
also play a crucial role to stabilize the size of the baryons. Indeed, without the CS term, we would 
be back to the scenario discussed in Ref. [7] in which the skyrmion, though stable and calculable, 
does not provide a good description of baryons. This CS term will at the same time be responsible 

^In Refs. [6,7] it was shown how the inclusion of the fuU tower of isovector resonances in a SU{2)l x SU{2)ii 
5D model can destabilize the skyrmion. 



1 



for the anomalous processes of the mesons such as tt — > 77, u — > 777, etc., which we will compute 
and compare with the experimental data. 

We claim that this is the first consistent holographic approach to baryons. Previous studies, 
although useful to understand the AdS/CFT correspondence for baryon physics, faced several 
problems. The first approaches considered only a truncated theory of resonances [11], and therefore 
had the same problems as 4D skyrmions. Later studies [6, 12, 13] were performed within the Sakai- 
Sugimoto model [14]. It was shown, however, that baryons are not calculable in this framework as 
their inverse size is of the order of the string scale which corresponds to the cut-off of the theory [6]. 
We will show that in our case we have an expansion parameter that, although not very small in 
real QCD {i.e. for = 3), allows for a pcrturbativc approach. Once a pcrtubative series has 
been formally built up, the accuracy of the predictions one gets at a given order depends on how 
fast the series converges, and not on the smallness of the naive expansion parameter. As we will 
see, the agreement with data of our leading-order predictions is quite good, suggesting that the 
series could converge fast in our case. Another important difference of our analysis is that we will 
find the non-linear soliton solution numerically, without the need of unreliable approximations. 



2 A five-dimensional model for QCD 

The 5D model that we will use to describe QCD with two massless flavors is the following. We 
will consider an U{2)l x U{2)r gauge theory with metric ds^ — a{zy {rjni.dx^dx'^ — dz'^), where 
represent the usual 4 coordinates and z, which runs in the interval [-^uv, -^ir], denotes the extra 
dimension. We will work in AdSs where the warp factor a{z) is 

aiz) = ^ , (1) 
z 

and zuy ^ to be taken at the end of the calculations. In this limit zi^ coincides with the AdS 
curvature and the conformal length 

fZlR 

dz. (2) 



We will denote respectively by Lm and Rm, where M = {fi,5}, the U{2)l and U{2)f> gauge 
connections. These are parametrized by = L\[aa/2 + Lm'^/2 and = -R^^o'a/2 + Rm'^/2, 
where are the Pauli matrices. The chiral symmetry breaking is imposed on the boundary at 
z = zi^ (IR-boundary) by the following boundary conditions: 

(L(U - Ryii) \z=ziYi. = ) (Lm5 + R^s) \z=zm. = ' (3) 

where the 5D field strength is defined as Tjmn = ^m'^n — d^'L'M — "^[Lm? Ljv], and analogously 
for Rmtv- At the other boundary, the UV-boundary, we impose generalized Dirichlet conditions 
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to all the fields: 



(4) 



We will eventually be interested in taking the 4D "sources" and to vanish. 



The AdS / CFT correspondence tells us how to interpret the above 5D model in terms of a 4D 
QCD-like theory, whose fields we will generically denote by \E'(a;) and its action by S^. This is 
a strongly coupled 4D theory that possesses an U{2)l x U{2)r global symmetry with associated 
Noether currents ^. If the 4D theory were precisely massless QCD with two flavors, the currents 
would be given by the usual quark bilinear, = Q''l,r1^Q\,r- Defining Z[/^,r^] as the 

generating functional of current correlators, the AdS/CFT correspondence states that [15, 16] 



PvE'exp 



I?LM2^RMexp[zS5[L,R]] 



(5) 

where the 5D partition function depends on the sources Z^, through the UV-boundary conditions 
in Eq. dl). 

The correspondence Eq. ([5]) leads to the following implication. Under local chiral transfor- 
mations, Z receives a contribution from the f/(2)^ anomaly, which is known in QCD. This 
implies [16-18] that the 5D action must contain a CS term 



CS 



'247r2 



j [uj,{L) - u,{R)] 



(6) 



whose variation under 5D local transformations which does not reduce to the identity at the UV 
exactly reproduces the anomaly. The CS coefficient, which is anyhow quantized even from a 
purely 5D point of view, will be fixed to A^c = 3 when matching QCD. The CS 5-form, defining 
A = -iAMdx'^, is 



cu5(A) = Tr 



A(rfA)2 + ^A^dA) + ^A' 
2 5 



(7) 



When A is the connection of an U (2) group, as in our case, one can use the fact that SU (2) is an 
anomaly-free group to write as 



o^5(A) = ^ATr [F^] + ^1 (dl)' + c^Tr 



AAF - -AA^ 
4 



(8) 



where A = A + Al/2 and A is the SU{2) connection. The total derivative part of the above 
equation can be dropped, since it only adds to Scs an UV-boundary term for the sources. 



^We are not considering the [/(l)-S'[/(A*'c)^ QCD anomaly, responsible for the 77' mass, since this is subleading 
in the large- A^c expansion. 
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The full 5D action will be given by S5 = Sg + Scs, where Sg is made of locally gauge invariant 
terms. Sg is also invariant under transformations which do not reduce to the identity at the UV- 
boundary, and for this reason it does not contribute to the anomalous variation of the partition 
function. Taking the operators of the lowest dimensionality, we have 

39 = -/ d'x j^^''^ dza{z) ^ |Tr [WL^ + y W^'^'^ + {L ^ . (9) 

We have imposed on the 5D theory invariance under the combined {x — x, L ^ i?}, where x 
denotes ordinary 3-space coordinates. This symmetry, under which Scs is also invariant, corre- 
sponds to the usual parity on the 4D side. Notice that we have normalized differently the kinetic 
term of the SU(2) and U{1) gauge bosons, since we do not have any symmetry reason to put 
them equal. In the large- A^c limit of QCD, however, the Zweig's rule leads to equal couplings (and 
masses) for the p and u mesons, implying a = 1 in our 5D model. Since this well-known feature of 
la.Tge-Nc QCD does not arise automatically in our 5D framework (as, for instance, the equality of 
the p and u masses does), we will keep a as a free parameter. The CS term, written in component 
notation, will be given by 

Scs = ^,j rf^x|ie^^^«^«LMTr[L;vo^PQ] + ^c^^^^^^LmL^vo^pq - {L ^ /?}| . 

(10) 

The 5D theory defined above has only 3 independent parameters: M5, L and a. 

From Eq. ^ we can extract the current operators through which the theory couples to the 
external EW bosons. These currents are obtained by varying Eq. ([5]) with respect to (exactly 
the same would be true for r^) and then taking Ip. = = 0. The variation of the l.h.s. of 
Eq. ([5]) simply gives the current correlator of the 4D theory, while in the r.h.s. this corresponds 
to a variation of the UV-boundary conditions. The effect of this latter can be calculated in the 
following way. We perform a field redefinition + 5L^ where 5L^(a;, z) is chosen to respect 

the IR-boundary conditions and fulfill 5L^(x, zuv) = "^^m- This redefinition removes the original 
variation of the UV-boundary conditions, but leads a new term in the 5D action, 55*5. One then 
has 

I j d^xTi [(j^(x))5/^(x)] = z j VUiVUm^S^ [L, R] exp [iS^ [L, R]] , (11) 

where the 5D path integral is now performed by taking l^j, = = 0, i.e. normal Dirichlet 
conditions. The explicit value of 5S^ is given by 

5^5 = J d'^xTr[Jl{x)Slp{x)] + j d^x{EOM) ■ 5L , (12) 

where Jl^ = -^E^cr" + Jifil and 

Jl^ = M,{a{z)Ll,) 1,=,^^ , J,^ = a'M,{a{z)L,^) . (13) 
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The last term of Eq. (fT2!) corresponds to the 5D "bulk" part of the variation, which leads to the 
equations of motion (EOM). Remembering that the EOM always have zero expectation value i we 
find that we can identify of Eq. ( fT3|) with the current operator on the 5D side: (j^)4D = (Ji,)5D- 
Notice that the CS term has not contributed to Eq. (fT2l) due to the fact that each term in Scs 
which contains a dz derivative (and therefore could lead to a UV-boundary term) also contains 
or R^ fields; these fields on the UV-boundary are the sources and that must be put to zero. 



3 Baryons as 5D skyrmions 
3.1 The static solution 

The QCD baryons correspond to the solitons of the above 5D theory, also referred as 5D skyrmions. 
These static solutions, exactly like the ones considered in Ref. [7], have unit topological charge 

B = / d^x / dze^ppaTr [L^'^L"" - R'^^R^"] , 

32^2 J J^^^ M-P- L J 

which we identify with the baryon number. The derivation of the solitonic configurations closely 
follows the one of Ref. [7]. There are, however, few differences. First, in Ref. [7] we were using 
time-reversal symmetry to put consistently to zero the temporal component of the gauge fields. 
Here, on the contrary, we have to use time-reversal t ^ —t combined with L —L and R —R, 
under which the CS is invariant. This transformation reduces, in static configurations, to a sign 
change of the temporal component of L and R and of the spatial components of L and R. We 
can therefore consistently put them to zero. For the remaining fields one can impose "cylindrical" 
symmetry, i.e. invariance under combined SU{2) gauge and 3D spacial rotations, and parity 
invariance under the combined action L ^ R and x — x. This determines our ansatz to be 



4"(x, z) = -Rti-^, z) , L»(x, z) = /?^(-x, z) , Lo(x, z) = /?o(-x, z) , (14) 



and 



_ 1 + 02(r,2:) Mr,z) ( 2. _ \ , 

So = (15) 

a r 

The above ansatz reduces our solitonic configuration to 5 real functions in 2D, Afi = {Ai,A2}, 
= 01 + ^02 and s, where x'^ = {r,z}. We notice that we have a residual U{1) invariance 



•^We have actually shown this here; notice that (5L^ was completely arbitrary in the bulk, but the variation of 
the functional integral can only depend on Slf^ ~ (5L^(a;, ^^uv)- 
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corresponding to = Qr 
charge +1 and s is neutral. 



exp[i a{r, z)x'^aa/ {2r)] under which is the gauge field, has 



The contribution of Sg to the energy is easily computed. The contribution from the SU{2) 
part is given in Eq. (11) of Ref. [7], while the U{1) part only adds the kinetic energy of s. We 
then have 



E„ = SttM, 



dr 



-21 R 



dz a(z) 



|2\2 



1 



(16) 



Notice that s has a negative kinetic term, since it corresponds to the temporal component of a 
gauge field. The CS term gives also a contribution to the energy. From the ansatz of Eqs. (|T^ 
and (flSj) we have LQ{r, z) = Ro{r, z) = s/ («r), that allows us to write the CS energy as a coupling 
of s to the topological charge density (the baryon number density), given in Eq. (12) of Ref. [7]: 



E, 



CS 



SttM. 



dr 



^IR 



where 
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dz-e""" 
r 



Nr. 



+ h.c.)+F: 



From Eqs. f[T6|) and ([T 



167r2 MgLa ' 
we can extract the EOM. One finds 



(17) 



;i8) 




'1 - 



{z){t<P*D,<P + h.c.)+^Le^'df,{^ 



(19) 



dj,{-i(t)*D^(t) + h.c.) + F-, 







The skyrmion configurations will be the solutions to these EOM with boundary conditions enforc- 
ing a definite topological charge B. For the B = 1 solution, these boundary conditions are given 
in Eqs. (13) and (14) of Ref. [7] for the SU{2) fields, while for the neutral scalar s we must impose 
s = at the three boundaries, z = zuy, r = and r oo, and impose d^s = at z = zjr. The 
IR and UV-boundary conditions come respectively from Eqs. (j3]) and (jl]); the one at r = arises 
from regularity, while the one at r ^ oo is necessary for the energy to be finite. The solution 
will be obtained numerically by the COMSOL 3.4 package [21]; a rescaling of the 0i^2 fields, as 
explained in Ref. [7], must be performed in order to avoid singularities at r = 0. 

At this point it is important to show that the 5D skyrmion configuration is stable and its 
properties can be consistently calculated within the 5D effective theory described above. This 
is easily established, along the lines of Refs. [6,7], in the case in which the CS term is a small 
perturbation to the gauge kinetic term, i.e., 7 <^ 1. In this limit, the solution is approximately 



''Any constant value for s would give finite energy, but since s = at the UV, taking s = is the only possibility. 
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given by a small 4D instanton configuration whose size p is determined by minimizing the energy. 
The AdS curvature induces a contribution to the instanton energy that grows as p/ L [7], while the 
CS term generates a Coulomb potential scaling as 7^i^^/ p^ [6]. Therefore the energy is minimized 
for p ~ By a Naive Dimensional Analysis (NDA) one gets @ A5 < 24:Tt^M^ and 7 ~ 

1/{K^L) <^ 1, so we have that p ^ I/A5 and therefore we expect the soliton to be insensitive 
to higher- dimensional operators. The situation is quite different in the Sakai-Sugimoto model [6]. 
There the energy dependence on p goes as -E(p) ~ p^/L^ + 7^L^/p^ where the first term comes 
from the curvature of the 5D space and the second from the CS with now 7 ~ l/(Mst-L)^, being 
Mst the string-scale. In this case one gets p ~ ^/iL ~ 1/Mst; the size of the soliton is of the 
order of the inverse cut-off of the model. This makes baryon physics totally sensitive to string 
corrections and therefore unpredictable, as already remarked in Ref. [6]. 

3.2 Soliton quantization 

Exactly like in the case of the 4D skyrmion [3] (see Ref. [4] for a comprehensive review), single- 
baryon states are described in our model as zero-mode time-dependent fluctuations around the 
classical static solution. Such zero modes, also called collective coordinates, are associated to the 
global symmetries of the theory; the situation is similar to that of the kink and the monopole 
[19]. The global symmetries of our static equations are U{2)v and 3-space rotations plus 3- 
space translations. The latter are associated with baryons moving with uniform velocity and 
can be simply ignored if one is only interested in static properties like magnetic moments, the 
axial charges and charge radii. The action of U{l)v is trivial on all our fields and then we 
are left with SU{2)y and 3-space rotations. The two rotations, however, have the same effect 
on our cylindrically symmetric solution so that rotating in 3-space would not lead to any new 
configuration which cannot be reached with only SU (2) rotation. Therefore, as in the case of the 
4D skyrmion, we only need to consider 3 collective coordinates which are encoded in a SU{2) 
matrix U. 

The zero-modes fluctuations we are interested in are constructed as follows. Let us perform 
an SU{2)v transformation U on the static solution discussed in the previous section. We obtain 

/?^(x, z;U) =U i?^(x, z) [/t , ^o(x, z; U) = i?o(x, z) , (20) 

^ In this NDA we are not considering the Nc dependence of 7. When this is included and our 5D model is 
matched to large-A^c QCD, we will see later that one gets 7^1 and p ^ L. Therefore, including the Nc factors 
will make the size of the baryon even larger with respect to the cut-ofF. 

^These quantities are defined in processes with very low transfer momentum during which the baryon only 
suffers a negligible acceleration. Our formalism, however, is non-relativistic so that uniform baryon motion cannot 
be ignored in a generic reference frame. The Breit frame is the correct one to compute form factors [4]. 
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where fi = 1,2,3,5; similarly for the U{2)l gauge fields. For a constant U, Eq. ( !20|) is also a 
solution of the EOM. We now introduce a small t-dependence on U, i.e. a small rotational velocity 
K* = — |Tr[cr*f/^(if//(it]. Eq. fl20|) is now an infinitesimal deformation of the static solution. In 
order for this configuration to describe a zero-mode fiuctuation, it should fulfill the time-dependent 
EOM to linear order in K and with c/K/rft = 0. Indeed, zero-modes correspond to directions in 
field space in which uniform and slow motion is permitted. 

We therefore need to solve the time-dependent EOM. Due to the invariance under the trans- 
formation {L i?, X <-> — x}, we can restrict the configurations to Li{'x,z,t) = —Ri{—:K,z,t), 
L5 o(x, z, t) = .Rs o(— X, z, t) and analogously for L, R. The EOM for the R fields, after separating 
temporal from spatial coordinates, read 



{a{z)R%) + ^-^e'^P^Rp^Rpa = 

Tr i^Ri/CjRpcT , 



~\~ 2 -^vCj R, 







a{z)DoR^^ 



jaL_ jli>pcr 



Do {a{z)R''^) 

ado (a{z)R'^'\ - aa{z)dQRQ^ - -fLe 



(21) 



flupcr 



Tr {RooRpa) + \RvoRpa 







where Euclidean metric is used to rise the spatial indices. We immediately see that, once the 
time-dependent ansatz in Eq. fl20l) has been chosen for the fields i?^ and i?0; the other components 
Rq and Rfi cannot be put to zero as in the static case. The time-dependence of U in Eq. (l20ll 
acts as a source for the latter components, as can be seen by looking at the first and the fourth 
EOM. As a result, Rq and Rj^ must be turned on. The same situation occurs in the case of the 
4D skyrmion of Ref. [4] in which the temporal and spatial components of the p and uj mesons 
are turned on in the skyrmion quantization (also in the case of the magnetic monopole for the 
temporal component of the gauge field). 

One can show that the second and the third EOM of Eq. (pT]) are solved, to linear order in 
K, by the ansatz in Eq. fl20|) if the fields Rq and Rf^ are chosen to be linear in K. Under this 
assumption, Roq is equal to BoRq up to terms proportional to dK/dt or quadratic in K that can 
be ignored in the approximation of slow motion discussed above. Furthermore, Roo and Rpa are 
of order K and DqR^^q = 0. Using this, the second and the third EOM of Eq. fl2T]) reduce to the 
static equations of sec. 13.11 



We are left with the first and fourth equations of Eq. (I2T1) . which are 7 elliptic equations for 
the 7 fields Rq and i?^; those are the analog of the Gauss law constraint for dyons in the case of 
monopoles [19]. To solve such equations we can again make a 2D ansatz following a generalization 
of the cylindrical symmetry we used in the static case in which the rotational velocity K also 
rotates together with the 3-coordinates x. The resulting 2D equations can be solved numerically, 
but we leave this for future work. 
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The quantization of the collective coordinates, from this point on, exactly proceeds as for the 
4D skyrmion [3,4]. First, one plugs the zero-mode configuration into the action and obtains a 
Lagrangian 

C = -M + ATr [d^UdQU^] , 

where M = Eg + Ecs is the classical mass of the soliton obtained from the static solution, while 
A depends also on the solutions for i?^ and Rq. Now, the collective coordinates U are treated 
as quantum mechanical variables, reducing the problem to the one of quantizing a spherical rigid 
rotor with momentum of inertia A. The quantization is therefore performed in a standard way 
and the energy eigenstates are interpreted as baryon states. One finds, as expected in the large- iVc 
limit, an infinite tower of baryons with increasing spin/isospin. More precisely, baryons are in the 
(p/2,p/2) representation of the spin/isospin SU{2) x SU{2) group; the first two levels p = 1,2 
are interpreted, respectively, as the nucleons and the A multiplet. We do not need to repeat 
this procedure in detail here, but only give the following useful relation that can be derived from 
quantization rules [3]: 

Tr [Ua^U^a""] = -- T , (22) 
3 

where S*" and 1°^ are respectively the spin and isospin operators. 



3.3 Static properties of baryons 



There are several static baryon observables which are independent of Rq and i?^ and therefore can 
be computed by knowing only the static solution of sec. 13. II together with the anstaz in Eq. ( |20l) . 
By looking at Eq. f|T3|) we see that the spatial components of the vector and axial currents, and 
the temporal component of the scalar one are independent of Rq and i?^ up to order or dK./dt 
terms. By plugging Eq. (l20l) into these currents one finds 



J. 



At 



J^ 



Vi 



J^ 



- M5 a(2;uv) 

- M5 a{zuY) 
-2 M5a(zuv) 



i k-^k ' 



F - 



zvv 



a 



d,,s 



2UV 



(23) 



The axial current is defined as Ja^ = Jr/i — Jifi, the vector as Jy = Jr^ + Jlij. and the scalar 
vector current (1/2 of the baryon number current) is Jy fi = 1/3 (^Jr^ + Jifij- From Eq. fl23|) . 
using Eq. (!22|) . we can extract the baryon axial coupling qa, the isovector magnetic moment /xy. 
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and the isoscalar electric charge radius g: 

N 1 1"^ 
97ra7 L Jq 



2UV 



97r«7 

4,5 = -^^/ rfrr^Kz)^..],^^^. (24) 



1 f°° 

— J drr^[a{z)D, 



7r7 



For the definition of these observables we follow the conventions used in Ref. [4]. It should 
be noticed that Eq. (!24|) has the right scaling in [20] if, in accordance with the AdS/CFT 
expectation, the 5D parameters 7, a and L scale like NJ^. We have numerically calculated all the 
quantities of Eq. fl2^ together with the mass of the baryon which is given by the classical mass 
M. The 3 parameters of our 5D model, M5, L and a, are determined from (1) the pion decay 
constant = 8M5 J dz/a{z) [10], (2) the mass of the p, nip ~ 3tt/{AL) [8,9] and (3) the decay 
constant ratio F^^/Fp = a. The results are given in Table [1] where we compare them with the 
experimental values. We restricted our comparison to nucleon observables (including the mass, 
even though all baryons are degenerate at the classical level), since we expect higher spin/isospin 
states (like the A-multiplet) to receive larger corrections. H We recall that these predictions are 
obtained in the semiclassical (large- iVc) limit, and therefore they should be valid up to corrections 
~ l/Nc- Consistently with this picture. Table [1] shows agreement with experiments within 30%. 





Experiment 


AdSs 


Mn (MeV) 


940 


1140 




0.79 


0.94 


9A 


1.25 


1.0 


fJ-p — fJ-n = 2//y 


4.7 


3.9 



Table 1: Predictions for the baryon static quantities where M5, L and a have been determined 
from the experimental values of F^ = 87 MeV, rUp = 775 MeV and F^/Fp = 0.88. 



In any model with exact spontaneously broken chiral symmetry the axial coupling is related 
to the pion- nucleon coupling Qt^nn by the Goldberger-Treiman relation 

This relation holds in our 5D model and can be used to derive Qt^nn from the value of qa and M^v 
obtained above. We get gnNN — 13.1 that is quite close to the experimental value 5'7rArAf|exp — 13.5. 



^In our model iS" — XK"" /2. Since we are performing a small K expansion, small spin (and isospin) states are 
predicted more reliably. 
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4 Meson anomalous couplings from the CS term 



The CS term of Eq. ([TO) is responsible for the anomalous couplings of the mesons to the photon 
and among themselves. □ In order to compute amplitudes involving a (real or virtual) photon at 
leading order in the electric charge e, one could proceed "holographically" and use directly Eq. ([5]) 
and Eq. (fT3l) to compute matrix elements of the electromagnetic current J^^ = + Jy^. 

Looking at the explicit form of this current, one easily realizes that electromagnetic interactions 
only proceed, in this "holographic basis", through the single exchange of resonances. Our model 
implements the Vector Meson Dominance (VMD) hypothesis, in which the photon only couples 
to hadrons by the mixing with the p and the uj. 

A second way to perform the same calculations, which in some cases is simpler, is to use 
a Kaluza-Klein (KK) expansion of the 5D gauge fields, in which the photon mixing matrix is 
automatically diagonalized. This can be done by giving a dynamics to the photon, which is the 
source associated to the electromagnetic current in Eq. ([5]). Making the source dynamical means 
integrating also over it in the path integral and this is the same as changing from Dirichlet to 
Neumann the UV-boundary condition of the corresponding 5D field + V^, where V = L + R. 
The electric charge is fixed to its experimental value e by adding a localized kinetic term for the 
photon of the form — (1/e^ — c.t.)F^/4 where the "counterterms" c.t. are needed to cancel the 
charge renormalization induced by QCD {i.e. bulk) effects. The latter are finite as long as ^uv is 
finite. At leading order in the charge e, however, the localized kinetic term is infinite (as e — 0), 
but the KK decomposition in the presence of a localized kinetic term is well known. One has a 
massless zero-mode with fiat wave function for V"^/3 + V^, that we identify with the "diagonalized" 
photon ^4^7'', which is now a mass eigenstate. One also has massive KK's which obey, in the limit 
e ^ 0, Dirichlet conditions at the UV. One basically returns to the original theory with Dirichlet 
5D fields, but with an extra photon zero-mode A|7'' which decouples as e — 0. The other relevant 
states are the pions vr'*, which are identified with the zero modes of the fifth component of the axial 
gauge bosons, and the p and uj resonances which are respectively the first isosinglet and isotriplet 
vector KK-states. We then have 

Al{x,z) = -±=f^{z)n%x) + ..., 

U^,^) = ^A^;^ + ^=fv{z)u;, + ..., 
o ^ ay/M^L 

V^{x,z) = ^A^y) + -^fy{z)p, + ..., (26) 
where fn{z) and fv{z) are respectively the 5D wave-function of the pions and vectors. We have 



*For a study of these processes in the Sakai-Sugimoto model see Ref. [23]. 
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Experiment AdSc 



T{uj TTj) 0.75 0.86 

T{lo Svr) 7.6 6.1 

r(p ^ n-f) 0.068 0.072 

r(u; ^ vr/i/i) 8.2 • lO""^ 7.9 • 10"^ 

r(io Tree) 6.5 • lO'^ 7.8 • lO^^ 



Table 2: Prediction of the anomalous partial decay widths in MeV where M5, L and a have been 
determined from the experimental values of = 87 MeV, nip = 775 MeV and F^/Fp = 0.88. 



U{z) = l/{a{z)N^) and, for the AdS metric, fy = zJi^Myz) / {NyL) with My ~ 37r/(4L) ~ = 
where in both cases A^^^ = / dz a{z)ff{z) / L. Inserting Eq. fl26|) into the CS term, we obtain 
the couplings between the pions and two vectors: 



vv 



ASn^F^ ^''^ 167r2 F^ 



pniT 



X 



167r2 F^ 
where F^"' = e'^^P'' Fp„/2, 

1 

9 pin 



a 



pu ' 



dz a{z) f^fv and 



a 



pu 



X 



dz fnf^ 



V 



(27) 



(28) 



The value of x turns to be very close to 1; for AdS we find x ~ 1.18. We will understand later 
why this is the case. We also want to remark that Eq. fl27j) shows an interesting relation between 
the wyiT (and P7vr) coupling and gpmr, the coupling of the p to two pions. This relation is fulfilled 
for any five-dimensional space. 



From Eq. fl271) we can calculate several meson partial decay widths. The first term of Eq. ([2] 
leads to the decay of the 7r° to two photons in accordance with the anomaly prediction. The decay 
widths r(ti; — > 717) and T{p n^y) arise from the second term of Eq. fl27|) . while r(ti; Svr) 
proceeds through virtual rhos. This latter process is dominated by the lowest 

state, the p, whose upTr coupling is given by the third term of Eq. ([2 



A[uj^{p) ^n%qo) + TT+{q+) + TT (g_)] 



Nc 9 



pTVTV 



X 



-pupa 



+ D{{q+ + qof)+D{{q^ + q,f 



(29) 



where D{p'^) = rnj^/ {rnj^ —p^)- The predictions for these partial decay widths are given in Table [21 
showing a very good agreement with the experimental data. 

The CS term also contributes to different pion form factors. For calculating form factors, 
however, it is more suitable to work in the holographic basis. We have seen that in this basis the 
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model exhibits the property of VMD. For example, the decay uj ^ ix'y proceeds as — * vrp^"^ — 717 
and similarly for rho decays. As in VMD models [24] , this allows us to derive the following sum-rule 
that relates the w^t: coupling, Qui-y,,, and the cupvr coupling, Qujptt'- 

^.7^ = E '"' (30) 
^ nipi^) nip 

where are the rho's decay constants that can be found in Ref. [8,9]. It is easy to verify 
Eq. ( 1301) . Using Eq. (p7|l . Eq. (1301) gives xgp.„T,Fp ~ nip that, since x ~ 1, implies gpT^Fp ~ mp. 
This latter equation was derived in Ref. [9] valid for any five- dimensional model. We can now 
easily calculate form factors. The 7r°77* form factor, F7r^(g^), at low Euclidian momentum is 
dominated by the exchange of the uj. We then have, as in VMD, F^^(g^) = m'^/{m1 + g^), where 
we have normalized the form factor as F^^^ylO) = 1. We obtain 

2 

F;^0= — , a=^~0.03, 31 

in perfect agreement with the experimental value [22]: a|exp — 0.032 ±0.004. This form factor was 
previously studied in Ref. [25]. We can also calculate the partial decay width V{ijj -n^jj,{ee)) 
that proceeds through a virtual photon. This process is proportional to the C(J7r7* form factor that 
in our model is simply given by F^-j^lp"^) = A{uj TT'~f)D{p'^), where A{uj TC'y) is the on-shell 
u 7T'~f amplitude. The prediction for this partial decay width is given in Table [2j 

Let us finalize this section with the following comment. In five-dimensional models arising 
from string theory, the effective gauge theory consists of the DBI term and the CS term. From 
the DBI term arises not only the kinetic term of the gauge bosons, but also higher-dimensional 
operators suppressed by the string scale. Since the anomalous couplings discussed in this section 
can only arise from the CS and not from the DBI, these couplings will not receive corrections from 
higher-dimensional operators. 

5 Global fit to mesonic and baryonic observables 

In the previous sections we have computed several properties of baryons and mesons, and we have 
shown that they agree with the experimental values within 30%. Nevertheless, in order to gain 
a better insight on the statistical significance of this approach, we will perform in this section 
a combined fit to many (baryonic and mesonic) physical observables. Our list of observables is 
presented in Tabled we have taken the physical quantities calculated in this article, together with 
other mesonic observables calculated in Refs. [8-10,18]. Our global fit is carried out by minimizing 
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Experiment AdSs Deviation 



r{uj Svr) 

T{p TTj) 

T{u) — > Tree) 



9A 




nil 



m 



■p 



0.068 
8.2 • 10-^ 
6.5 • 10-3 



6.9- 10-3 
-5.5 • 10-3 
940 



775 
1230 
782 
153 
0.88 
87 
6.0 



0.79 



1.25 
4.7 
0.75 
7.5 



6.2 • 10-3 
-6.2 • 10-3 
1180 



0.072 
7.4 • 10-^ 
7.4 • 10-3 



850 
1390 
850 
175 
0.90 
91 
5.4 



0.87 



0.98 
3.7 
0.82 
7.1 



+10% 
+13% 
+9% 
+14% 
+2% 
+5% 
-10% 
-10% 
+12% 
+25% 
+21% 

-21% 
-22% 
+10% 
-6% 
+5% 
-9% 
+14% 



Table 3: Global fit of mesonic and baryonic physical quantities. Masses, decay constants and 
widths are given in MeV. The RMS error of the fit is 16%. Physical masses have been used in the 
kinematic factors of the partial decay widths. 



the root mean square (RMS) error defined as 



where O^j.^ denotes the predictions of our model for any observable, Og^p its experimental value 
and A^p the number of predictions minus the number of parameters. In our case Np = 15. 

The statistical meaning of this procedure is the following. All our predictions are obtained 
at leading order in the (IGvr^MsL)-^ ~ l/Nc expansion, so that quite large quantum corrections 
are expected, which translate in quite large theoretical errors. Let us assume that all the observ- 
ables have the same relative theoretical error ^, i.e. AOp^.^ = ^Ogxp! looking at Eq. fl52]) one 
immediately sees that RMSE(M5,L, a) = ^a/xV^j where is the usual chi-squared variable 
constructed with errors AO^^.^. The minimum of Eq. fl32|) therefore represents the minimal value 
of the relative corrections ^ for which our model would successfully pass the test, i.e. for which 

= Np. All the observables included in the fit have an experimental error < 10%, which we can 
neglect because it is smaller, as we will see, than the final RMS error of the fit. 

The global fit gives 1/L ~ 350 MeV, M5L ~ 0.017 and a ~ 0.9 (7 ~ 1.27). As expected, the 




(32) 
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value of a from the fit is close to 1, as predicted in the large- Ai'c limit. Also it is worth noticing 
that the value of M5L comes to be quite close to the prediction arising from matching the current- 
current vector correlator to QCD at large momentum [8,9], M5L = Ai"c/(24'7r^) ~ 0.013. The RMS 
error of the fit is 16%. From the fit one can see that mesonic quantities are better predicted than 
the baryonic ones. Indeed, mesonic quantities alone give a fit with a RMS error around 10%, of 
order of the experimental errors. 

6 Conclusions 

Using a simple holographic model for QCD, we have been able to compute the properties of the 
baryons. We have seen that the presence of the CS term, needed to reproduce the QCD anomaly, 
is crucial to stabilize the size of the baryon at around the GeV. The contribution of this CS term 
to the mass of the baryon is as important as the leading term, implying the need for a fully 
numerical analysis of the non-linear solitonic configuration, not carried out before in the literature. 

We have calculated the axial coupling, the vector magnetic moment and the isoscalar electric 
charge radius of the nucleus as a function of the 3 free parameters of the model. These predictions 
have shown a good agreement with data -see Table [TJ Our approach can be extended to calculate 
other baryonic physical quantities that we leave for future work. 

We have also seen that the CS term is responsible for anomalous processes involving an odd 
number of pions, some of which we have explicitly calculated (see Table |2]) and showed to have 
an excellent agreement with the experimental data. Finally, we have done a combined fit of the 
mesonic and baryonic quantities predicted by this holographic model (Table [3l) that have shown 
an agreement with experiments of 16%. 

Some final comments on the important issue of calculability, i.e. on the dependence of our 
results on higher-dimensional operators. These operators are suppressed by the cut-off of the 
theory A5. Naive dimensional arguments say that the maximal value of A5 is determined by 
the scale at which loops are of order of tree-level effects. Computing loop corrections to the 
operator which arise from the term itself, one gets A5 ~ 247r^M5. Nevertheless, one gets a lower 
value for A5 from the CS term due to the Nc dependence of 7. Indeed, at the one-loop level, the CS 
term gives a contribution of order M5 to the F^ operator for A5 ~ 2Atx^ M^/ Nc^^ . Even though the 
cut-off scale lowers due to the presence of the CS term, we can still have, in the large- A^^c limit, a 5D 
weakly coupled theory where higher-dimensional operators are suppressed. For M5 ~ Nd (167r^L), 
as expected from the AdS/CFT correspondence, one has in fact A5L ~ nI^'^ —>■ 00. The same 
factor protects the physics of the 5D skyrmions which have size L. 
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